Unusual beats of the perpendicular-current giant magneto-resistance 
and thermopower of magnetic trilayers 
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Oscillations of the giant magnetoresistance (GMR) and thermo-electric power (TEP) vs. both 
the thickness of the non-magnetic spacer and also that of the ferromagnetic slabs are studied in the 
current-perpendicular-to-plane (CPP) geometry of magnetic trilayer systems, in terms of a single- 
band tight-binding model without impurities. The spin-dependent conductance has been calculated 
from the Kubo formula by means of a recursion Green's function method and the semi-infinite 
ideal-lead wires trick. Additionally the TEP is obtained directly from Mott's formula. In general, 
the thickness oscillations of the GMR and the TEP may have just one or two (short and long) 
oscillations. The long period, related to spectacular beats, is apparently of non-RKKY type. The 
TEP oscillations are strongly enhanced with respect to those of the GMR, have the same periods, 
but different phases and a negative bias. 



The phenomenon of giant magnetoresistance (GMR) 
in magnetic multilayers has been intensively studied both 
experimentally [|l]j|] and theoretically [||-|7j for more than 
five years now. After the paper was published it has 
become clear that one can expect a large magnetore- 
sistance even in systems which have no impurities and 
no structural defects. As pointed out in Q|, the GMR- 
oscillations have not only RKKY- type components, but 
may additionally reveal large particular non- RKKY os- 
cillations (see below). For the case of semi- infinite ferro- 
magnetic slabs, a period of those extra oscillations with 
the spacer thickness in the CPP (current perpendicular- 
to-plane) geometry has been shown in |Q| to originate 
from values of the in-plane wavenumber k\\ at which at 
least one spectral density vanishes at the ferromagnetic 
interface. 

In the present letter we study the CPP-GMR (current- 
perpendicular-to-plane) behaviour of layered systems of 
the type W/F\/ 'S/ 'F2/W ', where W stands for a semi- 
infite ideal lead wire, F\ and F2 for different ferromag- 
nets, generally with different thicknesses, and S for the 
non-magnetic spacer. Thus, apart from the leads, we are 
dealing with the usual magnetic trilayer systems. Our 
aim is to show that very pronounced beats may also, 
under some circumstances, appear as a function of the 
thicknesses of the ferromagnetic slabs. Besides the gi- 
ant magneto-resistance effect, we also calculate the cor- 
responding effect for the thermo-electric power ( TEP) . 

The calculation technique, we have developed, is based 
on the Green's function recursion method i it is close 
to that of H with one essential exception, namely instead 
of working with finite systems in the in-plane space (x, 
y directions), we have reduced the whole problem to one 
dimension by performing the Fourier transform in the in- 
finite x-y plane. Hence, our Green's functions fulfil the 



following equation: 
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where for the simple cubic lattice t z , z > = t5 z i tZ ±\ and 

A a z = V a (z) + 2t ■ [cos(k x ■ a) + cos(k y ■ a)} . (2) 

In eqn.(||), Ep is the Fermi energy; t (< 0) is the 
nearest-neighbour hopping integral; a is the lattice con- 
stant, and V a (z) is the spin-dependent atomic poten- 
tial, which we assume to be constant for given mate- 
rial. Hereafter, both \t\ and the lattice constant are taken 
as energy- and length units, respectively, which formally 
corresponds to t = —1 and a=l. The conductance of 
electrons of spin a is given by the well-known Kubo for- 
mula, which in turn may be expressed in terms of Green's 
functions. In our case the relevant Green functions are 
defined by eqn. (Q), and the conductance reads: 

8e 2 N x N v 
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Here the index i, which counts the z-planes, is arbi- 
trary, since the same amount of current flows through 
every cross-section. G" stands for the imaginary part 
of G, and the integration is over the 2-dimensional Bril- 
louin zone. We refer the reader to Refs. J|,^] for the 
detailed description how the Green functions of interest 
can be found from recursion equations, and we only men- 
tion that due to the ideal-leads trick, the first step of the 
iteration is analogous to what is known as the square-root 
terminator in the continuous fractions method (9). The 
prefactor in Eq. (^) contains the conductance quantum 
e /h as well as the cross-section area N x N y to which the 
number of active conductivity channels is proportional. 
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Due to this effective reduction of the problem to one 
dimension, we could easily compute, by a rapid recur- 
sion procedure, the Greens functions in eq. (0) for a large 
number of fc||-vectors and afterwards perform the integra- 
tion of cqn. (Q) very accurately by means of the highly ef- 
ficient special k-points method of Cunningham fllCfl . The 
GMR has been defined as 



GMR = 
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where the superscripts of Y refer to the parallel and 
antiparallel magnetization configurations of the two fer- 
romagnets, whereas the subscripts refer to the carriers 
considered. 

Before we performed original calculations, we have 
tested our programs by considering systems studied by 
Asano et al. in |j| , namely: (i) a non-magnetic multilayer 
/A/B/A/B consisting of 4 sublayers, each three mono- 
layers thick with alternating ±1 potentials, and (ii) their 
model for the FeCr superlattice: We have found very 
good agreement concerning both the CPP-GMR per unit 
area as well as the densities of states at the Fermi level, 
although the system of || has got only 12x12 atoms in 
the x-y plane and uses free boundary conditions. 

Our results for the CPP-GMR of the W/F 1 /S/F 2 /W 
systems are presented in the figures below. The plots cor- 
respond to the atomic potentials equal zero for both up- 
and down-spin electrons in the spacer and in the infinite 
lead wires. In the ferromagnets, the potentials for elec- 
trons with spin a = ±1 take the following values: Vi a and 
V 2a for the f ^-configuration and Via, and V 2 (_ a ) for the 
| j-configuration, where V\ a resp. V 2a and V 2 ^ a ) concern 
the slabs Fx and F 2 , respectively. (The actual values are 
given below and are also denoted in the figure captions). 
Additionally we have assumed a perfect matching of the 
minority bands of the ferromagnets with the spacer bands 
by putting K14. — V 2 [ = (cf. Q). Finally, unless oth- 
erwise stated, we have chosen the Fermi energy Ep=2.5 
well above the potential barriers, see below. 

It can be seen from Fig.l and Fig. 2 that the CPP- 
GMR oscillates as a function of the thickness rif of one 
or both ferromagnetic sandwiches with a short period of 
w 2 monolayers (ML). Additionally, in Fig.l, but not in 
Fig. 2, pronounced beats are observed with typical repeti- 
tion lengths of roughly 10 ML. Apparently those beats 
always appear, if the depth of the potential well corre- 
sponding to the ferromagnetic sandwich(es) with vary- 
ing thickness is Vaf = —1.8 (Fig.l) or —1 (not shown), 
whereas in case of V 2 ^ = — 2 the beats are absent (Fig. 2). 

In case of Fig. 3, where again V 2 ^ = —2, but where only 
the spacer thickness n s is varied, one observes an inter- 
mediate period of ~ 4.5 ML, whereas the beats do not 
appear. 
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FIG.l: CPP-GMR of the system 3F 1 /6S/n f F 2 between 
two infinite lead wires. Kit = —2 and Kst = — 1-8 in the ferro- 
magnets (all other potentials are 0). Ep=2.5. Similar results 
apply for the system nfF\/QS/njFi with Kit = V^f = —1.8. 
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FIG. 2 The same as Fig.l, but with Kit = v ^ = ~ 2 - 
Again similar results apply for the system njFi/6S/nfF\ 
with Kit — Vi\ = —2. 
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FIG. 3: The same as Fig.2, but with n/=3 on both sides, 
as a function of the spacer thickness n s . 
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In any case it can be easily seen that the GMR os- 
cillations decay always roughly ~ n ~f X ^ analogously to 
the findings of Mathon et al. concerning the oscilla- 
tions vs. the spacer thickness observed in [J7j . It results 
from these observations that oscillations and beats of the 
GMR, similarly as those of the exchange coupling, see 
p2| , are not specific for the non-magnetic spacer, but 
concern the ferromagnetic [Bslabs, too. 

In the following we try to identify the RKKY- and non- 
RKKY periods visible in our figures, using arguments 
along the lines of j7|JTl|] . 

In Mathon et al. observed roughly analogous be- 
haviour of the GMR as a function of the spacer thickness 
in a simpler system, namely a spacer imbedded between 
two identical semi-infinite fcrromagnets. As in our case, 
the authors of assume that within the spacer the po- 
tential V(z) vanishes, but concerning the ferromagnets 
their assumptions are different, since they assume the 
same vanishing potential for electrons with majority spin, 
but a positive potential V for electrons with minority 
spin, so that in their case the spacer represents a poten- 
tial well for minority carriers. (In contrast, in our case we 
have a negative potential for majority-spin carriers in the 
ferromagnets, but vanishing potential for minority carri- 
ers, so that in our case the ferromagnets represent quan- 
tum wells for majority carriers, which are interrupted by 
a barrier represented by our spacer and bounded from 
outside by the ideal leads, which have the same potential 
as the spacer). It was shown in [jjj that (i) a contin- 
uum of wavenumbers k z contributes to the conductivity 
r, i.e. the dependence on the varying thickness L was V ~ 
J dk x dk y c(k x , k y ) ■ exp[2 i k z (k x , k v ,Ef) L], and (ii) that 
dominating contributions can arise a) from wavenum- 
bers, where k z (k x ,k y ;Ef) is stationary as a function 
of (k x ,k v ), and b) from particular cut-off wavenumbers 
(k x , k y ), where c{k x , k y ) vanishes abruptly. These cut-off 
wavenumbers correspond in our case to unbound major- 
ity carriers which in z-direction have just the minimal 
kinetic energy in the wires to evade to infinity. The con- 
tributions under a) and b) are the "RKKY' and " non- 
RKK F-contributions" mentioned already above. 

At this place one should note in any case that the " non- 
RKKY wavenumbers" under b) do not play a role in the 
thickness-dependence of the exchange coupling of the fer- 
romagnets across the spacer, see pj . 

In both cases a) and b), one can estimate the important 
wave numbers k z from the asymptotic equation (valid ac- 
tually for large thicknesses rif and constant potential) 



k z (k x , ky', Ef) 
Here, 
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(5) 



(k x , k y ) 



in case a), with Ep=2.b and the values 
= (±7r, ±7r), where cos k x + cos k y is extremal, we 

get k^ 1 ' = 0.722 for the carriers with V a = (i.e. for the 
down-spin and spacer carriers in case of aligned ferro- 

( 2) 

magnets), and k z — 1.721 for the up-spin carriers with 
y T = -1.8 and k z 2) = 1.823 for V T = -2. The corre- 



sponding wavelengths can then be calculated as in |Tl]] 

from the expression X(p,q) = n/(pkz^ + qk z 2 ^), which 
yields A(0, 1) = 1.825 (V = -1.8) and =1.723 (V = -2), 
i.e. ~ the short period of two monolayers (2 ML) seen 
in the figures. It is interesting that in Fig. 3, where 
only the spacer thickness is varied, only a longer wave- 
length of ~ 4.5 ML is visible, which would correspond 
to A(1,0) = 4.35. However, these particular " RKKY- 
wavenumbers" would apparently not produce the beats 
in Fig.l, which have a much larger repetition length of 
« VQML. 

In order to get the beats seen in Fig. 1 and in Fig. 4 
below, one needs apart from the wave characterized by 

(2) 

k z another wave with a wavenumber very close to it. 
As there seems no way to get such a wave from eqn. 
(5) with V\ = —1.8, we suggest it to be of non-KKKY 
character. This means, it should come from electrons at 
E = Ef with values of (k x , k y ) ^ (±7r, ±7r). More insight 
into the nature of non-KKKY oscillations can perhaps 
only be gained by taking into account localized and reso- 
nant states in addition to the extended ones used for the 
sketchy estimations of the RKKY wavenumbers above: 
In fact, although the localized states do not contribute 
directly to the conductance, they modify the density of 
states at Ef and thereby influence the conductivity. 

In addition to the conductance T(Ef), where Ef is the 
Fermi energy, we also calculated the thermopower S: 
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FIG. 4: The same as Fig.l, with additional results for the 
" GMTEF 1 (Giant Magneto- TEP effect, lower curve), where 
TEP means "thermo-electric power". Note the shift of the 
TEP results by +1.75. 

S(E f ) is related to T(E f ) by Mott's formula, Q: 



S« = -^^(d/dEJlogr^E) 
3 e 



(6) 



at E = Ef, where a denotes the spin of the carriers, T 
the Kelvin temperature, fc^ Boltzmann's constant and e 
the elementary charge. 

Analogously to the GMR, see eqn. ([|), we define 
the "Giant Magneto- T£P-effect" GMTEP, where TEP 
means thermo-electric power, by 
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GMTEP = 



S 



rr 



all 

tJ -j- 



15 J. 



(7) 



It can be seen in Fig. 4 that the oscillations of the 
GMTEP are even stronger than those of the GMR, al- 
though they have a significant and negative bias and a 
different phase than those of the GMR. 

Finally it should be noted that the amplitude of the 
oscillations observed in the figures 1 and 2 is typically of 
the order of 1% to 2% only: Thus the oscillations can only 
be detected by very accurate calculations as the present 
one, and only with interfaces of suffcient quality: In fact, 
we have "smeared" the results for GMR(i) as follows 

GMR(i) -> (x/2) • [GMR{i + 1) + GMR(i - 1)] 

+ (1 - x) ■ GMR(i) , (8) 

with x = 0.5 and x = 0.25, respectively, and obtained 
instead of Fig.l the "smeared" Fig. 5, where the oscil- 
lations arc no longer visible for for strong interdiffusion 
(a; = 0.5), whereas for x = 0.25, they are roughly reduced 
by a factor (1/2) and can still be recognized. 



V a representing the exchange splittings). Similar results, 
with even stronger beats, have also been obtained for 
the "Giant Magneto- T£P-Effect" , where " TEP' means 
the Thermo-Electric Power. The beats of the GMR have 
been shown to remain visible, if there is only moderate 
chemical interdiffusion in the boundary layers. Finally, 
it has been shown that the short periods of 2 monolay- 
ers visible in the results arise from the usual "RKKY"- 
type mechanism, whereas the beats arise from certain 
"non-RKKY" cut-off contributions, analogous to those 
suggested in |7j]. 
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FIG. 5: The same as Fig.l, but smeared in accordance with 
eqn. (Eh with x = 0.25 (dashed line) and x = 0.5 (solid line). 

In conclusion, within an s-band tight-binding model 
combined with the Green's function recursion method, 
we have reduced the calculation of the CPP-GMR effec- 
tively to one dimension. We have computed the CPP- 
GMR for pure systems composed of two different ferro- 
magnets separated by a non-magnetic spacer and sand- 
wiched between ideal infinite lead wires. Unusual beats 
of the GMR vs. the thickness of the ferromagnetic slabs 
have been observed and shown to depend strongly on 
the kind of ferromagnetic slabs (i.e. on the potentials 
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